Collective motion in a Hamiltonian dynamical system 
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Oscillation of macroscopic variables is discovered in a metastable state in the Hamiltonian dynam- 
ical system of mean field XY model, the duration of which is divergent with the system size. This 
long-lasting periodic or quasiperiodic collective motion appears through Hopf bifurcation, which is 
a typical route in low-dimensional dissipative dynamical systems. The origin of the oscillation is 
explained, with self-consistent analysis of the distribution function, as the emergence of self-excited 
"swings" through the mean-field. The universality of the phenomena is also discussed. 
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Dissipative systems often show periodic, quasiperiodic, 
and chaotic temporal behaviors at a macroscopic level, 
when set at far away from equilibrium. They are de- 
scribed as low-dimensional dynamics, the discovery of 
which has marked an epoch of nonlinear dynamics studies 
in physics. Recalling that the involved degrees of free- 
dom are large, such macroscopic dynamics are a result 
of collective motion out of high-dimensional microscopic 
dynamics. Thus the collective motion is an important 
issue to be studied from both viewpoints of dynamical 
systems with many degrees of freedom and of nonequi- 
librium statistical mechanics. Indeed it has been inten- 
sively and extensively studied for systems consisting of 
a large number of chaotic elements over one and a half 
decades 0, 0, S Q , although the underlying microscopic 
dynamics in these studies have been restricted to dissi- 
pative chaos. 

Does such collective motion exist in closed thermody- 
namic systems, or Hamiltonian dynamical systems with 
many degrees of freedom? According to the standard 
belief in statistical physics, it seems rather difficult. In 
equilibrium, by definition, it does not exist. Also if the 
system rapidly relaxes to equilibrium, there exists only 
decay process. Thus long-term residence at a nonequi- 
librium state is at least required for such collective mo- 
tion. Indeed, there have been some studies on long- 
lasting metastable state in the relaxation to equilibrium, 
when some Hamiltonian dynamical systems are set at 
certain initial conditions [a, IS 0- However, no clear 
demonstration of macroscopic low-dimensionality in such 
metastable state has been reported yet, as far as the au- 
thors know. 

Here we discover such collective oscillation of macro- 
scopic (thermodynamic) variables in a closed Hamilto- 
nian system with many degrees of freedom. The oscil- 
lation is sustained over a long time, and indeed the du- 
ration increases with the system size, suggesting the di- 
vergence in the thermodynamic limit. Furthermore, the 
oscillation appears in a similar way to the bifurcation in 
low-dimensional dissipative dynamical systems, suggest- 
ing the low-dimensionality of the collective motion. In 



the present Letter, we report the essence of the phenom- 
ena, while the details will be soon reported elsewhere || . 

We adopt the Hamiltonian dynamical system of mean 
field XY model, or globally coupled pendula 0, , 

N 2 -. N N 
i=l i=l j=l 

All the N pendula interact with each other through phase 
difference. Each pendulum has two types of motion; ro- 
tation at a higher and libration at a lower energy. The 
equilibrium state is exactly solvable and uniquely deter- 
mined only with the total energy density U = H/N. Note 
that the system shows continuous phase transition. 

We focus on the dynamics of the variance of momen- 
tum, T(t), and the magnitude of the mean field, M(t), 
respectively defined as, 

N N 

r (*) = ^£ft(*) 9 , M(t)e*>t f I^e^W. (2) 
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These are macroscopic variables; indeed, they are nothing 
but the temperature and the magnetization of the system 
in equilibrium. Since Q yields the constraint, 

2U = T(t) + l-M(t) 2 , (3) 

we mainly discuss the dynamics of M(t) in the following. 
Using @, the equations of motion are described as single 
pendula interacting with the mean field: 

6j = Pi , Pj = ~M (t) sin(0j - 0(t) ) . (4) 

To study relaxation phenomena far from equilibrium, 
we prepare nonequilibrium initial states. In particular, 
we take the initial condition of {0i\ and {pi} given as 
follows. For a given total energy density U, assign an 
initial magnetization Mo = M(0). First, the initial dis- 
tribution of {8i} is set as a Boltzmann distribution of the 
equilibrium state determined, irrespective of U, with the 
magnetization Mq and the temperature T eq (Mo), where 
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FIG. 1: A typical time series of M(t). U = 0.69 0, M = 1. 
The dotted line indicates the equilibrium value. The abscissa 
axis is log scale. 



FIG. 2: A time series of M(t) in the metastable state. N 
10 6 (solid) and 10 5 (broken) for U = 0.5 and Mo = 0.9. 



T eq {M) is the equation of state in equilibrium: 
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where Zq(Mq) is a normalization. Next, the distribution 
of {pj} is set as a Maxwell distribution determined with 
the temperature T (U, M ) = 2U - 1 + M 2 to fulfill ©: 



fJp;U,M Q ) = 
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exp 
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where Z p (U,Mq) is a normalization. Here 9i and are 
given independently of each other. To sum up, the initial 
distribution is given as, 

f (e,p; U, Mo) = fg(0; M )f°(p; U, M ). (7) 

The present nonequilibrium initial conditions are con- 
trolled systematically with (U, Mo). Here, if Mo is the 
equilibrium value for the given U, then /o is exactly the 
equilibrium distribution. In this sense, the present initial 
states are connected to equilibrium on the (U, Mo) plane. 

Now we study the relaxation to equilibrium. A typical 
time series of M(t) is shown in FIG.Q] At the beginning 
(t ^ 0), M(t) decays almost exponentially. However, 
M(i), then, does not simply reach equilibrium but stays 
at around a larger value than the equilibrium value over 
a quite long time (t ~ 10 4 ). After the long interval, M{t) 
departs from the plateau toward equilibrium logarithmi- 
cally. Finally, M(t) reaches the equilibrium value, and 
fluctuates around there. As shown, the duration time in 
the plateau increases with N, and indeed for N > 10 5 , 
it is almost impossible to reach the equilibrium within 
numerical simulation (say t < 10 6 ). The result suggests 
that the metastable state is sustained in a macroscopic 
time scale in the thermodynamic limit (N — > oo). 

Metastability in the present model has been inten- 
sively investigated for a decade, especially by taking a 



rectangular (so called waterbag) initial momentum dis- 
tribution 0, Q . The metastable states exist only in the 
region just below the critical energy of the phase transi- 
tion, and there M(t) and T{t) take smaller values than 
those in equilibrium. This leads to a branch of negative 
specific heat, which is discussed as a reflection of ther- 
modynamic metastability, i.e. the local minimum of the 
free energy potential. On the other hand, the present 
metastable state takes larger values of M(t) and T(t) 
than those in equilibrium, and exists over a broader re- 
gion across the critical point than the negative specific 
heat branch. Furthermore, since no local potential mini- 
mum exists around there, the present metastable state is 
not a thermodynamically but a dynamically stable state. 
Thus the present metastability is a novel state. 

We study the metastable state in more detail. The 
closeup of the time series of M(t) at the state shows a 
periodic oscillation (FIG.EJ. The oscillation is not due to 
the finiteness of N. In fact, the oscillation becomes more 
apparent with increasing N, in strong contrast with the 
fluctuation around equilibrium that reduces to zero. The 
corresponding power spectrum density (FIG. OH shows a 
large peak (/ w 0.166) with its harmonic component, 
which remains sharp with increasing N. This indicates 
that the periodic motion survives even in the thermo- 
dynamic limit. Although in much longer time scale the 
amplitude of the oscillation decreases little by little, the 
decline becomes more gradual with increasing N, which 
implies that the periodic motion lasts permanently in the 
thermodynamic limit. Thus we have discovered the col- 
lective periodic motion. 

Besides the periodic motion, the macroscopic variables 
take various temporal pattern depending on the total en- 
ergy. In higher energy region, M(t) shows more com- 
plex collective motion than just a periodic oscillation. 
Its power spectrum (the inset of FIG.IHJ indicates that it 
is a quasiperiodic motion on a T 2 torus. In lower energy 
region, on the other hand, the system rapidly relaxes to 
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FIG. 3: The power spectrum density of M(t) for U — 0.5 
and 0.75 (inset). N = 10 4 (red), 10 5 (green), and 10 6 (blue). 
M = 0.9. 



equilibrium, without staying in any apparent metastable 
state, and there M(t) is almost stationary. 

Strictly speaking, there exist some fluctuation around 
each of the stationary, periodic, and quasiperiodic mo- 
tion, which remain finite even with increasing N . In 
fact, in the periodic region, the locus of the delay coor- 
dinate vector of M(t) in the embedded phase space [H| 
is not completely one-dimensional but has a finite width 
around the periodic motion. Similarly, even in the seem- 
ingly stationary region, the fluctuation around the fixed 
point does not decrease with N. This residual fluctua- 
tion is regarded as a result of high-dimensional collective 
motion || , as is observed in dissipative systems Q . 

Except for the fluctuation, the temporal behavior of 
the macroscopic variable changes as (almost) stationary 
— > (almost) periodic — > (almost) quasiperiodic with in- 
creasing U. This is regarded as a "bifurcation" of the col- 
lective motion. Here we note the similarity to the typical 
bifurcation in low dimensional dissipative dynamical sys- 
tems, fixed point — > limit cycle — > torus, through Hopf bi- 
furcations. Thus the present bifurcation of the collective 
motion is basically described as that of low dimensional 
dynamical systems, in particular by Hopf bifurcations. 

We investigate the bifurcation in more detail. We cal- 
culate the amplitude of M(t) against U, as shown in 
FIG. 0] The amplitude begins to increase above some 
critical energy Ub, with the approximate dependence of 
(U — Ub) 1 ! 2 '■ This suggests that the bifurcation from sta- 
tionary to periodic oscillation is a Hopf type. We have 
also studied the change of the behavior of the macro- 
scopic variables against Mo, which is another controllable 
for the initial distribution, and confirmed again the Hopf- 
type bifurcation from stationary to periodic oscillation. 

Next, we study the origin of the collective motion, in 
particular the periodic motion, from the viewpoint of the 
microscopic dynamics. Since it is almost impossible to 
investigate the trajectories of such large degrees of free- 
dom, we study one-body distribution function, f(9,p,t). 
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FIG. 4: The amplitude of M(t) against U, in the vicinity 
of the bifurcation point. N — I0 5 , Mo = 0.9. 25 samples 
are calculated for each U. The inset is the log-log plot of 
the same amplitude against U — Ub, where Ub = 0.35 is the 
"bifurcation energy" obtained from fitting. The dotted lines 
indicate M(t) oc (U - U b ) 1/2 . 



A snapshot of f(9,p,t) for the collective periodic motion 
when the phase of M(t) is zero (FIG.JSJ) shows a pair of 
concentrated density at ((9 — cf))/2-K,p) rj typo's, ±1), be- 
sides the expected high density at the center. The pair of 
peaks rotates clockwise along the separatrix keeping the 
localization, without diffusing out. Since such localized 
rotation is absent in equilibrium, it is important to ask 
its origin and relation to the collective motion. 

To answer the question, we consider an ensemble of N 
independent parametrically driven pendula: 

0j=P 3 , Pj = -M ext (t) sin 9j, (8) 

where M ext {t) = g + hs'milt. The equations are al- 
most equal to J2J except in that M ext (t) is externally 
applied, while M(t) is self-consistently determined in 
®. Indeed, let M int (t) be the "mean field" given by 
(J2J), and if Mi nt {t) — M ext (t), then the system is ex- 
actly the original Hamiltonian system. The Poincare 
section of © on fit = 0( mod 2ir), an equivalent of 
the standard map [l2|, yields islands of tori just at 
(9/2%, p) w (=p0.3, ±1) that correspond to the 1:1 res- 
onance on Q. Moreover, even though we first distribute 
the elements homogeneously, the axisymmetry is eventu- 
ally broken down, while the point symmetry is preserved 
due to the requirement of the dynamics; for example, in 
a snapshot, the first and third (the second and forth) 
quadrant are dense (sparse). 

On the basis of the above analysis, we give a self- 
consistent explanation on the collective periodic motion. 
Once a considerable number of elements exist in the 1:1 
island, they make a periodic and finite contribution to 
M(t). Then the periodic M(t) drives the elements, which 
in turn makes the 1:1 island stable. Moreover, the peri- 
odic drive also localizes other elements by breaking the 
axisymmetry, which also stabilizes the periodicity. Thus 
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FIG. 5: A snapshot of the contour of the one-body distribu- 
tion function, f(0,p,t), corresponding to t = 1001.2 of FIG. 
N = 10 5 . The dash-dotted line indicates the separatrix. 



the periodic collective motion stably continues with the 
self-organization of self-excited pendula, or "swings." 

The above picture also explains the bifurcation from 
the stationary to the periodic state as follows. The 1:1 
resonance state requires quite a little energy. In low en- 
ergy region, only a small number of elements can have 
enough energy to stay in the 1:1 island, which does not 
contribute to the stability of M(t). In high energy region, 
on the other hand, considerable elements become able to 
have enough energy, which leads to stable periodic mo- 
tion of M(t). Thereby, with increasing U, the bifurcation 
to periodic motion appears. The bifurcation along the 
Mo-axis follows the same scenario; as Mq increases with 
constant U, Tq obtained from J2J also increases, whereby 
the elements have enough energy to stay in the 1:1 island. 

In summary, we have discovered collective periodic and 
quasiperiodic motion at a macroscopic level in the Hamil- 
tonian dynamical system of mean field XY model. In 
the thermodynamic limit, these motions are suggested 
to continue forever, as their lifetime seems to diverge. 
The collective motion appears through Hopf bifurca- 
tion, which is a typical route in low-dimensional dissipa- 
tive dynamical systems. Thus we have observed macro- 
scopic low-dimensional motion out of microscopic high- 
dimensional Hamiltonian chaotic dynamics. The mech- 
anism of the macroscopic periodic motion is explained 
from the viewpoint of the underlying microscopic dynam- 
ics, as the self-organization of self-excited swings. 

The above mechanism can work for Hamiltonian dy- 
namical systems with global coupling, or under mean 
field approximation, when the resulting one-body dynam- 
ics have the separatrix of motion, which allows for the 
self-consistent resonance of mean field oscillation. The 
existence of such separatrix in the one-body dynamics 
means the phase transition in terms of statistical physics. 
Hence it is expected that the collective motion we dis- 
cussed here can be universally observed in a system with 
relevant mean field interaction and phase transition. In 



fact, we have recently discovered that the mean field <fi 
model does show collective periodic motion as well as 
Hopf bifurcation, as in the present case @ ■ Since systems 
with global or long-ranged coupling with phase transition 
broadly exist in nature, for example in molecular clusters 
and (bio)polymer systems, it will be an important issue if 
the collective motions are observed in such real systems. 

In the end, we add a comment on Ref. 13] that have 
reported "coherent oscillating cluster" in the metastable 
state of the negative specific heat branch in the mean field 
</> 4 model. Although it may be a kind of collective mo- 
tion, it probably corresponds to the residual fluctuation 
in the almost stationary motion rather than the collective 
oscillation in our result, since the oscillation amplitude 
is quite small. Actually, in the present XY model, some 
tiny peaks in the power spectrum are also observed at 
the stationary motion when Mq = 1 8] . It is also noted 
that the Hopf bifurcation as in low-dimensional dynami- 
cal system is confirmed in our case. In fact, we have also 
discovered such low-dimensional collective motion even 
in the mean field (f> model, with initial distribution as in 
the present Letter rather than that adopted in Ref. [l3| . 
Finally, the present oscillation is observed, not in the neg- 
ative specific heat branch, as is also in contrast with the 
gravothermal oscillation 14] in a self-gravitating system. 

The authors are grateful to Kensuke Ikeda for discus- 
sion. This work was supported by a Grant-in-Aid for 
Scientific Research from MEXT Japan. 



* Electronic address: morita@complex.cu-tokyo.ac.jp 
Present address: Faculty of Science and Engineering, 
Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 
169-8555, Japan 

[1] K. Kaneko, Phys. Rev. Lett. 65, 1391 (1990), Physica D 
55, 368 (1992). 

[2] H. Chate and P. Manneville, Prog. Theor. Phys. 87, 1 
(1992). 

[3] G. Perez and H. A. Cerdeira, Phys. Rev. A 46, 7492 

(1992); N. Nakagawa and Y. Kuramoto, Prog. Theor. 

Phys. 89, 313 (1993); A. S. Pikovsky and J. Kurths, 

Phys. Rev. Lett. 72, 1644 (1994). 
[4] T. Shibata and K. Kaneko, Phys. Rev. Lett. 81, 4116 

(1998); T. Shibata, T. Chawanya, and K. Kaneko, Phys. 

Rev. Lett. 82, 4424 (1999). 
[5] M. Antoni and S. Ruffo, Phys. Rev. E 52, 2361 (1995). 
[6] See, for example, recent papers; V. Latora, and A. Rapis- 

arda, Physica D 193, 315 (2004); Y. Y. Yamaguchi et al, 

Physica A 377, 36 (2004); and references therein. 
[7] H. Morita and K. Kaneko, Europhys. Lett. 66, 198 

(2004). 

[8] H. Morita and K. Kaneko, in preparation. 
[9] T. Konishi and K. Kaneko, J. Phys. A 25, 6283 (1992). 
[10] S. Inagaki, Prog. Theor. Phys. 90, 577 (1993). 
[11] N. H. Packard et al., Phys. Rev. Lett. 45, 712 (1980); 
F. Takens, Detecting Strange Attractors in Turbulence, 
in Lecture Notes in Mathematics (Springer, New York, 



■5 



1981), vol. 898. [14] J. Makino, Gravothermal Oscillation, in Dynamical Evo- 

[12] B. V. Chirikov, Phys. Rep. 52, 263 (1979). lution of Star Clusters, (Kluwer, Dordrecht, 1996). 

[13] T. Dauxois, S. Lepri, and S. Ruffo, Nonlin. Sci. Num. 

Sim. 8, 375 (2003). 



